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Abstract. We analyze in details the statistical signif- 
icance of the claim by Bird [2002] of a power law dis- 
tribution of plate areas covering the Earth and confirm 
that the power law with exponent 0.25±0.05 is the most 
robust and parsimonious model for all plates, includ- 
ing the very largest plates, when taking into account 
the constraint that the plates areas must sum up to 4-7T 
steradians. We propose a general class of fragmentation 
models that rationalize this observation and discuss the 
implications for the earth dynamics and the general self- 
organization of tectonic deformations at multiple scales. 



Introduction 

According to the theory of plate tectonics, the out- 
ermost 100 km of the solid Earth is composed of a rel- 
atively small number of internally rigid plates (litho- 
sphere) that move over a weak substrate (asthenosphere) . 
The plate sizes and positions change over time and are 
driven by the internal mantle convection. DeMets et 
al. [1990] performed a global inversion to determine 
the relative rotation rates of the 12 largest plates (the 
NUVEL-1 model) later refined into the NUVEL-lA so- 
lution [DeMets et al., 1994]. The edges of these plates, 
where they move against each other, are sites of in- 
tense geologic activity, such as earthquakes, volcanoes, 
and mountain building. These plate edges are not sharp 
narrow boundaries but are often constituted of complex 
systems of competing faults and other geological struc- 
tures of width extending over several hundreds of kilo- 
meters for transform and subduction boundaries and up 
to thousands kilometers for continental collisions. It is 
now common lore to view such tectonic deformation as 
possessing some kind of self-similarity or fractal proper- 
ties [Scholz and Mandelbrot, 1989], or better a hierarchi- 
cal structure [Ouillon et al, 1996]. In addition, several 
researchers have repeatedly proposed to describe this 
multi-scale organization of faulting by a hierarchy of 
blocks of multiple sizes [Gelfand et al., 1976; Gorshkov 
et al., 2001]. In these models, blocks slide against each 
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other, rotate, lock at nodes or triple junctions which 
may represent the loci of major earthquakes, in a way 
similar to (but at a reduced scale compared to) the rel- 
ative motion of the major tectonic plates. An example 
of a complex network of faults is observed in the broad 
San Andreas transform boundary between the Pacific 
and the North American plates in California, for which 
Bird and Rosenstock [1984] have suggested 21 possible 
microplates within southern California alone. 

Keeping in mind these ingredients of a few major 
plates at large scales on one hand and a hierarchical 
self-similar organization of blocks at the boundary scale 
on the other hand, the recent reassessment of present 
plate boundaries on the Earth by Bird [2002] is par- 
ticularly inspiring: taking into account relative plate 
velocities from magnetic anomalies, moment tensor so- 
lutions, and/or geodesy, to the 14 large plates whose 
motion was described by the NUVEL-1A poles, model 
PB2001 [Bird, 2002] includes 28 additional small plates, 
for a total of 42 plates. Bird [2002] suggests that the 
cumulative- number/area distribution for his model fol- 
lows a power-law for plates of less than 1 steradian of 
area. 



Statistical analysis of the size 
distribution of the 42 plates 

Analysis with rank-order statistics and the 
Pareto law 

The statistical analysis of such a data set is very dif- 
ficult due to its small size N = 42 but is not impos- 
sible. Figure 1 shows the complementary cumulative 
number N(A) of plates as a function of area A in stera- 
dians, i.e., the number of plates with an area equal to 
or larger than A. Most of the data except for a few 
largest plates follow a linear dependence in the double 
log-scale of the figure. The slope of this straight line 
is close to 0.25. We first compare this sample with a 
pure power law (the Pareto distribution). For this pur- 
pose we use the rank-ordering method (see [Sornette et 
at, 1996]). We put the sample into descending order 
Ai > A2 > ... > An- The probability density function 
(PDF) of the n-th rank, denoted (f) n< jsr(x), is well-known 

4> n , N (x) = (N-n+l) (^j F N - n (x) (l-Fix))"- 1 f(x) 

(1) 

where F(x), f(x) are the distribution function (DF) and 
PDF of the random values in question. Putting the 
Pareto law F(x) = 1 — (a/x) 11 , x > a, into (|l|), we get 

^ N (x) oc (1 - {a/xYf- 1 x ^- n+ x)-i . (2) 



The mode M nt jy of the PDF (||) (i.e., the maximum 
of the PDF) is the most probable value of the random 
variable A n : 

M^.afSe^Y" . , 3) 
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Besides, an interval around the mode containing some 
prescribed probability (say, 90%) can be derived from 
the density 4>n,N{x). The dependence (||) with [i = 1/4 
is shown as the long-dashed line in Fig. 1. The two 
short-dashed lines represent the upper/lower limits of 
the 90%-confidence interval. The data are well ac- 
counted for by the power law, except, perhaps, for the 
three smallest ranks, i.e. the three largest plates, the 
Pacific, the Antarctica, and the Africa plates, which fall 
outside the confidence interval. 



"Finite-size" constraint on the PDF of plate 

sizes 

From a visual inspection of figure 1, it might be ar- 
gued that the deviation from the power law prediction 
(|J) occurs somewhat earlier, say, at rank n = 7, i.e., the 
seven largest plates with area more than 1 steradian be- 
long to a different population than the rest of the plates. 
The simplest possibility is that the largest plates are 
sensitive to a "finite-size" effect, namely that the sum 
of areas over all plates must sum up to An. The largest 
continent and ocean plates are commensurable with the 
size of the convection cells in the upper mantle and they 
would thus feel the constraints due to the finite size of 
the Earth surface. In contrast, one could argue that the 
smaller plates appeared as a result of the interaction 
and collision between the larger plates and do not feel 
any superficial restriction. More precisely, the smaller 
plates may be the result of fragmentation/creation pro- 
cesses in the neighborhood of the borders and especially 
near triple points of the seven largest plates. 

The distribution gN,c(x) of sample values xi, ...,xjv 
conditioned by the constraint 

Sn = x\ + ■■■ + %n — C , (4) 

where C is a constant (Att for the plates) is modi- 
fied from its unconditional Pareto density expression 
f(x) — [i a^/x 1+t - L for x > a and f(x) — for 
x < a. The lower threshold (the minimum value) 
for the plate data is a = 0.00178 steradian and cor- 
responds to the smallest documented plate in [Bird, 
2002] . Denoting the unconditional density of the sum 
Sk = x\ + ... + Xk, k = N — 1, N by Sk(x), we have 
9n,c(x) = s N ^x(C - x) f(x)/s N (C), for a < x < C. 



Thus, the constraint (Q) decreases the unconditional 
Pareto density f{x) by a factor sn-i(C — x)/sn{C) 
which acts as a "taper." 

In order to use the Maximum Likelihood (ML) method 
for the estimation of the exponent /i, we need the vec- 
torial distribution of conditional sample to take into 
account the interdependence between the different vari- 
ables (areas of the plates) Xi,...,X43 induced by the 
constraint ([|). The corresponding likelihood function 
is therefore 

. . S{xx + ... + X42 -4tt) , 
0(xi, ...,X42 MJ = 7T~\ J H1U -/142 W ■ 

(5) 

The resulting ML estimate is = 0.25 ±0.05. With this 
value, we generate an artificial sample of 1045 condi- 
tional 42-dimensional (42D) vectors with the condition 
(|J) with C — 4ir. Rank-ordering each of these 1045 vec- 
tors, we determine their sample medians Mi, M42, 
where Mj is the median of the j-th rank. These condi- 
tional medians are slightly smaller than given by for 
the unconditional Pareto distribution. The conditional 
distribution (|^) allows us to construct a confidence do- 
main for the 42D random vectors, defined as a "corri- 
dor" of the form [cMj] (l/c)Mj\, the constant c = 0.244 
being chosen such that 95% of vectors fall within this 
corridor. The medians Mi,...,M42 and their corridor 
are shown in Fig. 2: all samples of the tectonic plates 
falls within the 95% confidence corridor, showing that 
(in contrast with the "pure" Pareto used in Fig.l) the 
Pareto model together with the total area constraint 
accounts satisfactorily for all the data, including the 
largest plates. 

We have performed many other tests which all con- 
firm the goodness of fit of the constrained Pareto distri- 
bution (or of its approximation by the truncated Pareto 
distribution). Using Wilks' theorem [Rao, 1965], we 
used the so-called Gamma-distribution with comple- 
mentary cumulative distribution function G(x), trun- 
cated from both sides G(x\a, b) = C(a, &)§t+f for uq < 
x < u\ and find that (1) the introduction of the ad- 
ditional parameter a in the exponential does not in- 
crease significantly (probability of rejection 38%) the 
likelihood as compared with the truncated Pareto and 
(2) the truncated Pareto distribution is much preferable 
compared with the truncated exponential distribution 
(significance 1 — 10~ 8 ). This conclusion is strongly con- 
firmed by tests using so-called sufficient statistics for 
both the power law and exponential models. We also 
tested the popular lognormal DF, which is a natural 
model for fragmentation since it is the limit distribu- 
tion for processes of subdivision of rock under rather 



general assumptions (basically multiplicative processes) 
[Kolmogorov, 1941]. Using the Kolmogorov test and the 
bookstrap method, we estimate that the data deviates 
significantly from the lognormal DF at the 95% signifi- 
cance level. 

Discussion 

Our main result is that Bird's suggestion [2002] of a 
power-like distribution of plates areas is confirmed up 
to the largest plates, when taking into account the very 
natural constraint that the plates areas must sum up to 
4w steradians. 

Given this strength of the validity of the Pareto dis- 
tribution for the 42 plates documented in [Bird, 2002], 
one can expect that it can be extrapolated beyond this 
range to smaller yet undetected plates. Using the com- 
plementary Pareto distribution {a/xY with a = 0.00178 
steradian and fi — 0.25, this extrapolation predicts 100 
plates larger than 1.35 10~ 4 steradian (5,500km 2 or 
74 x 74 km) [Bird, 2002]. The total area needed to 
define another 58 plates ranging from this size up to 
the smallest size in PB2001 would be only about 0.028 
steradians, which could be taken from large plates like 
EU and NA without materially affecting their areas. As 
discussed by Bird [2002], the places where additional 
small plates are most likely to be recognized is within 
the zones of distributed deformation identified in [Bird, 
2002], which have total area of 0.838 steradians (6.7% 
of Earth). 

The main suggestion inspired by Bird's discovery is 
that the separation of scales between a plate tectonics 
involving a few major plates at large scales and tec- 
tonic deformation at smaller scales may be an artifact 
that resulted historically from insufficient statistics and 
more recently from a not quite correct interpretation of 
"visual" deviations from the pure power law. With the 
improved extended compilation of [Bird, 2002] which 
will probably still improve in the years to come, we 
conjecture that there is a hierarchy of plates and block 
sizes from the largest ones to the smallest scales, which 
may be the result of the self-organization of the plate 
motions, creation and destruction over millions of years. 
In other words, we suggest that our view of plate tec- 
tonics should be revised in terms of a dynamical model 
of plates with creation, fragmentation and destruction 
that may occur at all scales. This view may provide 
a link between the large organization of major plates 
down to the block structures at scales of tens of kilome- 
ters or less discussed in the introduction. 

What could be the mechanism(s) at the origin of the 
observed power law distribution of plate areas? We pro- 



pose a fragmentation model with source terms, which 
can be progressively enriched by taking into account 
the kinematic compatibility constraints between plates. 
There is a large class of such fragmentation models 
that can be solved explicitely within the scaling the- 
ory of [Cheng and Redner, 1988]. In a nutshell, their 
approach is to write down a linear integro-differential 
coupled equation for the time evolution of fragment 
populations. The fragment distribution p(l) is char- 
acterized by the moments m a — J Q l a p(l)dl. Note 
that calculating a moment corresponds to taking the 
Mellin transform of p{l). The moments m a give ac- 
cess to p(l) by taken the inverse Mellin transform. The 
equations of evolution for the fragment population then 
transform into the following recurrence equation for the 
moments: m a+ p = lo l~" l m ai where j3 is the exponent 
describing the overall rate of breakup (assumed to be 
proportional to I 13 ), uj is a normalizing constant and 
L a = Jq x a b(x)dx, where b(x) is the probability that 
the fragment ratio be x. From the factor L _ 1 , if there 
is a value a* such that L a * = 1, then all moments with 
a > a* will become infinite. Provided reasonable an- 
alyticity conditions hold, it follows that the value a* 
is a pole of m a . Taking the inverse Mellin transform 
of m a then allows us to get p(l) and, using the exis- 
tence of the pole at a*, this immediately predict that 
p(l) ~ l^( 1 +l 3+a "> is a power law distribution. This 
model was used by Ouillon ct al. [1996] to account for 
the power law distribution of joint spacing in a granitic 
massif in Saudi Arabia, with exponent fj, s = 1/2. Note 
that block sizes are proportional to spacing taken to 
the second power, which then predicts fi = (J, s /2 = 1/4, 
which is in agreement with our best estimate /i = 0.25 
for the plate area distribution. This agreement be- 
tween two very different systems may be fortuitous but 
is worth keeping in mind, especially in view of our claim 
of a continuity of the physics of fragmention from the 
large plates to the blocks. We note also that the frag- 
mentation model of [Cheng and Redner, 1988] predicts 
generically log-periodic corrections to the main power 
law since, for a large variety of b(x), there are solutions 
of L a * = 1 with complex exponents a*. This has been 
shown to signal the existence of discrete scale invari- 
ance [Ouillon et al.. 1996]. It is tempting to interpret 
in this way the three approximately regular oscillations 
observed in figure 1, with a prefered scaling ratio ap- 
proximately equal to 10. This value is in good agree- 
ment with the square of the prefered scaling ratio 3.5 
found in the discrete hierarchy of linear sizes in fracture 
of rocks by Sadovskiy et al. [1984]. 

Our present conclusion on a universal fractal char- 



acter of plate fragmentation does not prevent the exis- 
tence of some underlying distinctions. For instance, P. 
Bird (private communication and [Bird, 2001]) proposes 
three main tectonic origins for plates: plume-assisted 
rifting for the larger plates (ranks 1-8 in Figs. 1 and 
2), continental collision stress for the intermediate plate 
sizes (ranks 8-20), and back-arc spreading for the small- 
est plates (ranks 21-42). In Figs. 1 and 2, one can dis- 
cern slightly different patterns of the tail behavior of 
these three subgroups. However, any formal statistical 
analysis of these possible distinctions would be, to our 
opinion, excessive due to the extreme smallness of the 
data set. 
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Figure 1. Complementary cumulative distribution of 
the areas of tectonic plates (open circles) compared 
to the fit with the formula (ft) for a power law (cen- 
tral long-dashed line) with exponent fi = 0.25 and 
a = 0.002. The small-dashed line and medium-dashed 
line provide a confidence interval defined in the text. 



10 



95% confidence domain 




Plate area, sterradian 



Figure 2. Medians Mi, M42 (continuous line) and 
their corridor at the 95% confidence level delimited by 
the two dotted lines of the conditional Pareto distri- 
bution cf>(xx, ...,X42\fJ') given by (Q) compared with the 
empirical cumulative distribution of the 42 plate areas 
(circles linked by straight segments), documented by 
Bird [2002]. 
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